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Work  during  the  grant  period  from  June  15,  1978  through 
June  14,  1979  included:  (1)  the  measurement  of  time  resolved 
electron  densities  in  laser  triggered  discharges,  (2)  the  study 
of  the  initial  buildup  of  charge  and  current  in  these  gaps,  and 
(3)  improvements  of  the  data  acquisition  equipment  used  for  this 
work.  An  effort  was  also  directed  towards  studying  spectroscopi¬ 
cally  the  events  occurring  between  the  time  of  arrival  of  the 
triggering  laser  pulse  and  the  formation  of  the  arc. 

1.  Experimental  Apparatus 

For  all  discharge  studies  the  experimental  arrangement  shown 
in  Fig.  1  was  used.  The  spark  gap  was  enclosed  in  a  vacuum  con¬ 
tainer  which  could  be  evacuated  and  then  backfilled  with  gas  to 
the  desired  pressure.  One  electrode  of  the  gap  had  a  small  hole 
drilled  in  it  to  allow  the  introduction  of  the  triggering  laser 
beam  into  the  gap  in  a  longitudinal  geometry.  For  all  experiments 
reported  here  the  electrodes  were  made  of  aluminum.  The  triggering 
laser  was  a  Molectron  UV-400  nitrogen  laser  which  delivered  5  mJ  • 

O 

pulses  of  3371  A  radiation  in  10  ns.  The  roughly  rectangular  (spa¬ 
tially)  output  beam  was  focused  to  a  point  on  the  lower  electrode 
with  a  quartz  lens,  through  a  CaF2  window  in  the  vacuum  cell. 

Emission  from  the  gap  passed  through  another  CaF2  window  and 
was  collected  and  focused  onto  the  entrance  slit  of  a  0.5  m  spec¬ 
trograph  by  a  quartz  lens.  After  passing  through  the  spectrograph, 
the  radiation  was  focused  onto  the  photocathode  of  an  SIT  vidicon 
detector  connected  to  an  optical  multichannel  analyzer.  This  setup 

O 

allows  a  spectrum  covering  a  range  of  200  A  to  be  acquired  in  a 


single  shot  with  good  sensitivity.  The  intensifier  section  of  the 
SIT  may  be  gated,  allowing  time-resolved  data  to  be  acquired  with 


Schematic  drawing  of  experimental  setup. 


a  temporal  resolution  cf  down  to  50  ns.  For  these  time- resolved 
studies,  the  gate  generator  could  be  triggered  after  an  adjustable 
delay,  either  from  the  nitrogen  laser  or  from  the  rising  edge  of 
the  gap  current. 

Voltage  was  applied  to  the  gap  through  a  charged  cable  arrange¬ 
ment.  A  length  of  RG-8  cable  was  charged  to  the  operating  voltage 
through  a  large  charging  resistor  connected  to  a  regulated  high  vol¬ 
tage  power  supply.  Application  of  the  triggering  pulse  to  the  gap 
then  resulted  in  a  nicely  shaped  current  pulse  of  duration  propor¬ 
tional  to  the  length  of  the  charged  cable  (about  1  vs  for  most  expe¬ 
riments)  .  In  order  to  minimize  waveform  distortion  due  to  reflec¬ 
tions,  the  gap  was  terminated  in  a  carefully  matched  50  n  load,  and 
the  gap  current  was  determined  by  measuring  the  voltage  across  this 
load  resistor. 

2 .  Initial  Buildup  of  Charge  and  Current 

A  typical  oscillogram  of  the  gap  current  in  a  laser-triggered 
discharge  in  H2  at  350  Torr  is  shown  in  Fig.  2.  The  observed  current 
rises  smoothly  to  a  maximum  value  corresponding  to  the  complete  col¬ 
lapse  of  the  gap  voltage.  This  behavior  is  in  contrast  to  that 
observed  under  somewhat  similar  conditions  with  conventionally  trig¬ 
gered  gaps,  where  well-defined  current  steps  or  even  maxima  are 
observed  prior  to  the  complete  closure  of  the  gap.  These  features 
have  been  associated  with  the  formation  of  the  cathode  fall,  the 
arc  channel  formation,  and  finally  the  bridging  of  the  cathode  fall 
region  by  the  arc  channel  [1-5].  Apparently,  in  laser  triggered 
gaps,  the  sequence  of  events  is  significantly  different  than  in 
more  conventionally  triggered  devices. 
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TIME  (200  ns/div.) 


CURRENT  PULSE  IN  LASER-TRIGGERED  GAP. 
H2  at  350  Torr,  5  kV,  50  ft  load. 

Figure  2. 
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The  primary  difference  between  the  two  triggering  methods 

is  related  to  the  amount  of  charge  liberated  by  the  trigger. 

Conventional  U.V.  trigger  sources  provide  a  comparatively  small 

2  3 

number  of  electrons  (10  -  10  )  whereas  laser  triggering  re¬ 

sults  in  a  much  larger  number.  In  order  to  characterize  better 
the  initial  pulse  of  positive  and  negative  charge  provided  by 
the  triggering  laser,  initial  experiments  were  conducted  in 
which  the  gap  was  evacuated  as  completely  as  possible  (<0.1 
mTorr)  and  a  small  voltage  in  the  range  of  0  to  20  V  applied 
to  the  gap.  Under  those  conditions,  charge  multiplication  from 
ionization  and  other  secondary  processes  should  be  minimized, 
and  the  integrated  gap  current  should  be  a  measure  of  the  quan¬ 
tity  of  free  charge  in  the  plasma  fireball. 

The  results  of  these  experiments  proved  quite  interesting. 
Two  current  waveforms  observed  under  these  conditions  are  shown 
in  Fig.  3.  Figures  3a  and  3b  show  the  current  observed  with 
plus  and  minus  23  V  across  the  gap  respectively.  For  the 
polarity  in  Fig.  3a,  the  positive  charge  is  drawn  across  the 
gap,  whereas  in  Fig.  3b,  the  negative  charge  traverses  it. 

In  spite  of  the  striking  qualitative  differences  in  the  cur¬ 
rent  pulse  shapes,  the  total  charge  in  each  case  is  nearly  the 
—  8  —8 

same  (2.8  x  10  and  2.4  x  10  C, respectively) ,  supporting 

the  contention  that  charge  multiplication  does  not  play  an 

important  role  in  these  results  and  that  the  total  charge  in 
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the  fireball  is  in  the  range  of  2  -  3  x  10  C. 
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This  quantity  of  charge  can  have  significant  effects.  For 
example  if  the  charge  were  unshielded  and  confined  to  a  sphere  of 
1  mm  radius,  the  electric  field  at  the  surface  of  the  sphere  would 
be  of  the  order  of  2  x  106  V/cm.  Of  course,  in  these  experiments?  . 
such  fields  are  not  produced,  because  the  plasma  remains  approximately 
neutral,  but  it  is  clear  that  space  charge  effects  play  an  important, 
and  probably  even  dominant  role  in  our  laser  triggered  experiments. 

In  light  of  the  important  role  space  charge  effects  play  in 
breakdown  phenomena,  a  theoretical  effort  to  understand  them  was 
recently  initiated  and  is  to  be  part  of  the  continuing  program  of 
research.  In  a  simple  model,  the  positive  and  negative  charge  carrier 
densities,  p+(r,t)  and  p_(r,t),  are  expected  to  obey  the  following 
relations  [6] 


3  P  +  ->■  |  ->■  I  2 

__  =  _  u  +  (  P+E)  -aU_p_  |e  |  +D+V  p+ 

3  P  .  2 

-r—  =  -y_v*  (  p_E)  +  ocy_p_  |e  |  +  D_  7  p_ 


(1) 


where  p+  and  p_  are  the  densities  of  positive  and  negative  charge, 
y+  and  D+  are  the  mobilities  and  diffusion  constants,  a  is  the  Town¬ 
send  ionization  coefficient  (a  function  of  if)  ,  and  £  is  the  electric 
field,  including  space  charge  effects.  The  charges  are  assumed  to 
drift  in  the  electric  field  and  to  diffuse  in  response  to  density 
gradients.  Positive  and  negative  charges  may  be  created  by  ioniza¬ 
tion  of  gas  molecules  (Townsend  a  processes) .  Charge  creation  at 
the  electrode  surfaces  appears  as  a  boundary  condition.  Although 
not  specifically  included  here,  other  effects  such  as  photoionization, 
electron  attachment ,  or  negative  ion  production  could  be  added  as 
necessary.  Since  the  electric  field  depends  on  the  charge  density, 
both  the  drift  term  and  the  Townsend  a  term  make  these  equations 
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non-linear,  and  very  difficult  to  solve  analytically. 

Efforts  toward  developing  a  computer  program  to  solve  these 

I 

equations  numerically  are  continuing.  In  this  work,  symmetry  about 
the  axis  of  the  discharge  is  assumed,  but  radial  variations  are 

specifically  included.  Kohrmann  and  others  [6-9]  have  reported  the 

I 

results  of  somewhat  similar  calculations  in  which  the  discharge  was 
assumed  to  be  uniform  in  the  radial  direction,  thus  reducing  the 
problem  to  a  one-dimensional  one.  However,  such  an  assumption  is 

<: 

clearly  not  justified  for  laser-triggered  discharges,  and  radial 
non-uniformity  is  essential  for  proper  description  of  the  phenomena. 
Yoshida  and  Tagashir^and  later  Davies,  Evans,  Townsend,  and  Woodison^ 
have  reported  results  of  two  dimensional  calculations.  Because  of 
the  two-dimensional  nature  of  the  problem,  considerable  care  must 
be  taken  to  minimize  computer  execution  time,  and,  particularly, 
memory  requirements.  The  problem  is  most  acute  in  the  calculation 
of  the  electric  field  from  the  charge  densities.  At  the  end  of  the 
grant  period,  we  had  not  found  a  satisfactory  solution  to  the  pro¬ 
blem  of  calculating  the  electric  field  which  required  reasonable 
computer  resources  while  providing  acceptable  accuracy.  Work  is 
continuing,  however. 

3.  Electron  Density  Measurements 

Electron  density  information  is  important  for  understanding 
electrical  breakdown  processes.  Accordingly,  Stark  broadening  of 
the  Ha  and  Hg  lines  was  utilized  to  obtain  time  dependent  electron 
density  information  from  laser-triggered  discharges  in  hydrogen. 

In  addition  to  measuring  the  electron  density  during  the  discharge, 
we  were  also  interested  in  determining  densities  after  the  termination 
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of  the  arc.  For  this  purpose,  a  voltage  source  which  could  be 
rapidly  turned  off  at  a  set  time  was  needed.  The  charged  coaxial 
cable  system  described  earlier  functioned  well  in  this  regard.  By 
careful  adjustment  of  the  load  resistor  a  clean  current  pulse,  as 
shown  in  Fig.  2,  with  no  ref lection- induced  after-pulses  was  obtained. 

Using  a  50  ns  wide  gate,  time  resolved  spectra  of  the  Hg  line 
were  obtained  at  a  number  of  times  during  and  after  the  current  pulse. 
In  the  initial  experiments,  the  gating  generator  was  triggered  by 
the  laser.  Because  of  jitter  in  the  delay  between  the  triggering 
laser  pulse  and  the  leading  edge  of  the  gap  current  pulse,  it  was 
necessary  to  record  a  spectrum  and  a  current  pulse  oscillogram  simul¬ 
taneously  on  a  single  shot  basis.  In  order  to  obtain  a  spectrum 
at  a  specific  time  in  the  current  pulse,  it  was  necessary  to  monitor 
several  such  shots  until  one  was  obtained  with  a  current  pulse  pro¬ 
perly  lined  up  with  the  gate  pulse. 

Besides  making  data  taking  tedious,  this  problem  made  it  diffi¬ 
cult  to  signal-average  over  a  number  of  shots,  and  it  made  it  essen¬ 
tially  impossible  to  obtain  the  data  necessary  for  using  Abel  inver¬ 
sion  techniques  [12-14]  to  unfold  the  spatial  variation  of  the  elec¬ 
tron  density.  In  order  to  overcome  this  difficulty,  the  gating 
generator  was  triggered  from  the  leading  edge  of  the  current  pulse. 
This  procedure  essentially  eliminated  the  temporal  jitter  but,  due 
to  an  internal  delay  of  about  100  ns  in  the  gate  generator,  measure¬ 
ments  during  the  early  parts  of  the  leading  edge  of  the  current 
pulse  could  not  be  obtained. 

Figures  4  and  5  show  typical  spectra  of  the  H  line  seen  at 

P 

various  times  in  the  current  pulse  shown  in  Fig.  2.  The  first 
spectrum  was  obtained  from  a  single  shot,  the  second  is  the  sum  of 
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20  shots.  Variation  of  the  linewidth  is  quite  evident.  Following 
the  analysis  of  Griem  [12],  the  electron  density  is  approximately 
given  by 

Ne  =  C(Ne,T)yx  3/2  (2) 

c  s 

where  C  is  a  constant  which  is  only  weakly  dependent  on  the  plasma 
temperature.  A  plot  of  the  density  thus  obtained  is  shown  in  Fig. 6. 

It  should  be  noted  that  the  data  in  Fig. 6  represent  spatially 
averaged  densities.  Because  of  the  geometry  of  the  collection  optics 
system,  high  spatial  resolution  is  available  in  the  direction  paral¬ 
lel  to  the  plane  of  the  entrance  slit  of  the  spectrometer,  but  there 
is  essentially  no  depth  resolution  so  that  the  spectra  obtained 
represent  data  averaged  over  this  dimension  of  the  discharge.  Assum¬ 
ing  that  the  discharge  has  cylindrical  symmetry,  the  radial  variation 
of  the  spectra  may  be  unfolded  by  stepping  the  arc  image  across  the 
entrance  slit  of  the  spectrograph  and  obtaining  a  spectrum  for  each 
step.  With  these  data,  an  Abel  inversion  may  be  performed  to  obtain 
the  spatially  resolved  spectra  [12-15].  For  this  purpose,  an  effi¬ 
cient  computer  program  compatible  with  the  limited  capabilities  of 
our  Tektronix  4051  minicomputer  was  developed  for  performing  numeri¬ 
cally  the  Abel  inversion  of  experimental  raw  data.  A  brief  descrip¬ 
tion  of  the  program  follows. 

Experimentally,  the  image  of  the  arc  is  moved  across  the  entrance 
slit  of  the  spectrograph  and  a  spectrum  is  obtained  for  each  of  a 
regularly- spaced  net  of  image  positions.  If  is  the  number  of  fre¬ 
quency  points,  or  channels,  in  a  single  spectrum,  and  Nx  is  the  num¬ 
ber  of  spectra  acquired  in  sweeping  the  arc  image  across  the  entrance 
slit,  then  this  procedure  produces  an  Nx  x  data  matrix  which  cons¬ 
titutes  the  input  to  the  inversion  program,  In  n  . 


If  for  a  given  frequency,  the  functional  form  of  the  lateral 
variation  of  the  intensity,  I(x),  is  known,  then  the  radial  varia¬ 
tion  of  the  intensity,  e(r)  is 


e(r) 


1_  I'(x) 

*  *4  /x^  _  r2 


d  x 


(3) 


where  r0  is  the  radius  of  the  arc.  Because  of  the  difficulties 
associated  with  numerical  differentiations,  for  each  frequency, 
the  N  data  points  representing  I(x)  were  least-sauares  fit  to  a 
sum  of  Chebyshev  polynomials,  yielding  an  approximation  to  the 
lateral  intensity  function, 


y~  a jTj (x)  =  I (x)  (4) 

j 

The  necessary  derivative  and  integral  operations  in  eq  (3)  could 
then  be  carried  out  analytically,  and  the  radial  emission,  e(r) 
evaluated  for  a  set  of  Nr  radial  positions,  ra. 

The  coefficients,  a..  in  eq  (4)  are  most  easily  evaluated  if 
the  input  data  are  spaced  with  a  cosinusoidal  as  opposed  to  a 
regular  distribution.  Quadratic  interpolation  was  used  to  con¬ 
struct  such  a  data  set  from  the  regularly-spaced  input  set.  With 

these  considerations,  the  radial  emission  intensity  matrix  becomes 

Nx  Nx  L 

S*i  1  ji 

ven 


where : 


1^  is  the  input  intensity  for  frequency  at  lateral  point 

X3. 

Skj  is  an  interpolation  matrix  for  converting  the  regularly 
spaced  data  set  to  a  cosinusoidal  spacing. 

Tah  *  Tl*xk^  where  Tn  is  the  Ath  Chebyshev  polynomial. 
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and  L  is  the  highest  order  included  in  the  polynomial  sum.  Only 
even  terms  are  included  in  the  l  summation  to  ensure  symmetry 
about  the  discharge  axis.  can  be  computed  by  writing  T4  as 

an  explicit  sum,  with  the  result6 
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where  k  =  A/2  -  j  -  1. 

For  all  wavelengths,!^,  the  matrices  S,  £,  and  B  remain  constant 
and  therefore  need  be  calculated  only  once.  Therefore  if  C  is 
defined  by 

**  *  Si  i  b«  (7> 

%  even 

then  the  radial  emission  matrix  is  simply 


'1,  -  Hi  Cl,  V 


(8) 


Evaluation  of  this  expression  for  all  spectral  points  (500  in 
our  case)  and  Nr  radial  points  (typically  about  20)  is  then  rela¬ 
tively  fast,  but  the  volume  of  data  which  must  be  manipulated  severely 
taxes  the  limits  of  the  4051,  and  it  is  necessary  to  use  the  magnetic 
cassette  tape  drive  for  intermediate  data  storage  during  the  cal¬ 
culation. 
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The  inversion  algorithm  was  tested  by  analytically  constructing 
a  synthetic  lateral  intensity  function  I (x)  from  a  chosen  functional 
form  for  e(r).  This  synthetic  I (x)  was  fed  into  the  inversion  rou¬ 
tine  and  the  output  e(r)  compared  with  the  original.  For  all  cases 
tested  the  agreement  was  excellent. 

Figure  7  and  8  show  the  Abel  inverted  spectra  for  a  range  of 
radial  positions  in  the  arc.  Fig. 7  results  from  data  taken  near 
the  start  of  the  arc  phase  (t=150  nsec.),  and  Fig. 8  from  data  taken 
just  after  the  current  shut-off  (t=1.2  usee.).  The  integrated 
spectral  intensity  of  the  Hc  line  is  plotted  in  Figs  9  and  10  for 

P 

these  same  two  times.  These  curves  were  obtained  by  drawing  straight 
lines  through  the  experimental  points,  and  therefore  they  accurately 
represent  the  noise  in  the  signal  acquisition  and  manipulation  pro¬ 
cesses.  From  these  data  it  is  clear  that  the  luminous  arc  channel 
is  expanding  with  an  average  velocity  of  approximately  5xl04  cm/sec. 

The  electron  density,  Ne,  was  found  from  the  experimental  Hg 
spectra  using  eq  (2) .  The  electron  densities  obtained  in  this 
manner  for  two  different  times  are  shown  in  Figs  11  and  12.  The  tem¬ 
perature  assumed  for  these  data  was  20,000°K,  but  the  results  would 
not  have  been  greatly  affected  if  a  temperature  of  up  to  40,000°K 
had  been  assumed.  In  these  results,  the  edge  of  the  arc  channel 
is  clearly  visible  at  a  position  consistent  with  the  edge  of  the 
luminous  column  as  seen  in  Figs  9  and  10.  Additionally  a  second 
peak  in  the  electron  density  is  seen  further  out.  We  believe  this 
second  peak  results  from  a  shock  wave  travelling  radially  outward 
from  the  channel  center.  Shock  waves  have  been  observed  very 
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clearly  in  experiments  using  Schlieren  photographic  techniques. 
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but  these  results  mark  the  first  clear  observation  of  a  shock  wave 
in  the  electron  density.  More  work  needs  to  be  done  to  put  the 
shock  wave  interpretation  on  a  firmer  foundation  and  to  better 
characterize  it. 

4.  Early-Time  Fluorescence  Measurements 

An  effort  was  made  to  study  spectroscopically  the  early-time 
sequence  of  events  starting  with  the  small  plasma  fireball  produced 
by  the  triggering  laser.  Unfortunately,  in  hydrogen,  very  little 
emission  from  the  gap  is  observed  during  these  times  [1]  so  that 
these  efforts  have  been  unsuccessful  to  date.  Several  emission 
lines  associated  with  the  electrode  material  were  seen,  but  no 
emission  which  could  be  identified  with  either  atomic  or  molecular 
hydrogen  was  observed. 

Much  of  the  molecular  hydrogen  emission  spectrum  is  in  a  con¬ 
tinuum.  Measurement  of  the  emission  spectrum  with  a  lower  disper¬ 
sion  spectrograph  should  be  useful,  therefore,  in  enhancing  the 
detected  intensity.  A  coarser  grating  has  recently  been  acquired 
for  this  purpose.  Also,  other  gases  such  as  N2  are  reported  to 
have  a  brighter  emission  during  the  early  phases  of  breakdown,  and 
they  should  be  investigated  also. 

5 .  A  Microprocessor-Based  Interface 

The  PAR  optical  multi-channel  analyzer  which  we  are  using  is 
designed  to  be  used  in  conjunction  with  a  small  computer  which 
serves  to  control  it  and  to  store  and  manipulate  data.  A  Tektronix 
4051  minicomputer  along  with  a  locally  constructed  interface  unit 
are  currently  being  used  for  this  purpose.  Since  the  4051  is  rela¬ 
tively  slow,  during  a  run  the  interface  must  accumulate  and  store 


t; 


data  from  the  analyzer  for  later  transfer  to  the  4051.  Although 
the  present  interface  unit  performs  adequately  in  this  function, 
greater  flexibility  is  needed  both  to  allow  simple  arithmetic  ope¬ 
rations  such  a  background  subtraction  to  be  performed  within  the 
interface,  and  to  allow  for  memory  expansion  and  for  more  efficient 
use  of  the  memory.  Better  use  of  memory  is  especially  important 
for  applications  requiring  two-dimensional  scanning  of  the  OMA 
vidicon. 

Accordingly,  a  microprocessor-based  interface  unit  was  designed 
and  is  being  built.  The  unit  uses  the  Z-80  microprocessor  and 
offers  the  following  features: 

(1)  Memory  may  be  expanded  by  simply  plugging  in  more 
memory  cards  and  modifying  the  Z-80  program. 

(2)  Data  may  be  stored  in  words  of  one  to  three  eight-bit- 
byte  lengths. 

(3)  Data  transfers  to  the  4051  are  facilitated  by  putting 
the  data  in  the  correct  internal  format  of  the  4051. 

(4)  The  contents  of  memory  may  be  displayed  in  graphic 
form  on  an  oscilloscope  or  chart  recorder  without 
going  through  the  4051. 

(5)  Arithmetic  operations  on  data  stored  in  memory  may 
be  accomplished  internally. 
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